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For certain systems, the N-particle ground- state wavefunctions of the bulk happen to be exactly 
equal to the N-point space-time correlation functions at the edge, in the infrared limit. We show 
why this had to be so for a class of topological superconductors, beginning with the p+ip state in 
D=2+l. Varying the chemical potential as a function of Euclidean time between weak and strong 
pairing states is shown to extract the wavefunction. Then a Euclidean rotation that exchanges 
time and space and approximate Lorentz invariance lead to the edge connection. We illustrate 
straightforward extension to other dimensions (eg. "^He- B phase in D=3+l) and to correlated 
states like fractionalized topological superconductors. 

PACS numbers: 



The boundaries or edges of condensed matter sys- 
tems received scant attention until recent developments 
showed them to be fertile areas of research both in the 
Fractional Quantum Hall Effect (FQHEjP^. and in topo- 
logical insulators and superconductors-^. 

In two spatial dimensions, the edge dynamics is de- 
scribed by conformal field theor}^' w hich w as also used 
to produce wave functions in the bullff^^. Moore and 
Read^ showed that one may view the FQHE wavefunc- 
tions and the quasi-hole excitations as conformal blocks 
in which both electrons and the quasiparticle coordinates 
are treated on the same footing and their charges and 
braiding properties are severely constrained. For an ex- 
haustive review of many related topics see Nayak et a¥^. 

What are the minimal ingredients necessary to estab- 
lish equality of edge correlations and bulk wavefunctions? 
Are analytic functions or d=2 conformal invariance re- 
quired? We show that our edge-bulk equality follows for 
a class of topological superconductors in various dimen- 
sions invoking only approximate Lorentz symmetry. The 
connections obtained here using an effective low energy 
hamiltonian differ from CS theorjOni [yi which the hamil- 
tonian vanishes and only non-dynamical particles enter 
via Wilson loops, as reviewed in Ref, 13, 

We shall first write down an operator expression for 
Z(J), the generating function of A^-body wavefunctions 
of the bulk. This is shown to be accomplished by intro- 
ducing a time dependent chemical potential that changes 
abruptly at some Euclidean time. We then drop some 
high derivative terms which do not matter in the infrared, 
and express Z{J) as a Grassmann integral over a Lorentz 
invariant action. Rotating by 90 degrees to exchange 
time and a spatial direction we obtain the same topo- 
logical superconductor but with a spatial edge induced 
by the jump in chemical potential. We find that the 
same Z{J) has now morphed into the generating func- 
tion for the correlation functions of the edge excitations. 
Three examples are given: the p -\- ip superconductor in 
I) = 2 + 1, ^ He B phase in D = 3 -\- 1 and a p-wave 
superconductor (the Ising model) in D = 1 + 1. 

Extracting Wavefunctions: Recall that given a 
second-quantized A'-body state |$) with wavefunction 



(/)(xi,X2, ..xn) we extract <p using 

(I){xuX2,..Xn) = (0|^(xi)...^(xAr)|<l>). (1) 

where (0| is the Fock vacuum and ^ is the canonical elec- 
tron destruction operator. For problems with variable 
number of particles, let us define the generating function 

Z{J) = (0|e/^^'^(^)^(^)|<l>) (2) 

which yields A'-body wavefunctions upon differentiating 
A"- times with respect to the Grassmann source J{x). 

We want to express Z{J) as a path integral when |<l>) is 
the ground state of a Hamiltonian H without conserved 
particle number. Since Euclidean time evolution for long 
times projects to the ground state, we can obtain |$) as 

m = uio-,-^)\i) (3) 

where \i) is a generic initial state and /7(0~, — oo) is the 
imaginary time propagator from — oo to 0~. Then we in- 
sert the operator exp [J J{x)'^{x)dx] at time 0. Finally, 
we obtain the Fock vacuum by evolving a generic state 
(/I from time +oo to 0+ using a hamiltonian H' with 
a huge negative fi that empties out fermions so that we 
may write (0| = (/|/7((X), 0+). Thus 

Z{J) = (/|/7(oo,0+)e^^(")^(")^"/7(0-,-oo)|z) (4) 

which has a path integral representation. 

Example 1: p-\-ip : The mean-field hamiltonian i^I^E^ 

Tr_\^fJ^ ^( ak^-fi A-{ki-ik2)\f Ck 

(5) 

here 1, 2 are spatial indices and xs will be time. We em- 
ploy the minimum k dependence in the pairing function, 
and set the coefficient A = 1 for convenience so the gap 
function is: A(/ci, k2) = ki — ik2 

The ak'^ term is needed to ensure the nontrivial topol- 
ogy of the weak-coupling phas^^^ and to populate it with 
electrons for /i > 0. We shall remember this association 
but drop the '/c^' term in the computations since it does 
not affect infrared correlations. 
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FIG. 1: (a) Wavefunction: The original superconductor with 
= /i- > Ues in the xi — X2 plane and evolves in Euclidean 
time X3 from — oo to 0~, projecting out the ground state |$). 
At xs = 0^ the chemical potential drops abruptly to a large 
negative value leading to the Fock vacuum, (b) Corre- 
lation functions: A Lorentz rotation makes xi the new time 
and X3 a the spatial coordinate along which the system has 
an edge at xs = 0. The world-sheet of the edge lies in the 
xi — X2 plane at xs = 0. 



Now the mean field Hamiltonian in real space: 

H=jd?x ^^(-/i)^ + ^(^^(-^^1 - ^2)^^ + h.c) 

^(6) 

leads to corresponding Grassmann action for [/(O, — oo): 

/oo 
d?x / dx^[i)ViP ^ {l)idi) ^ ^IjiB^p] (7) 
-co J — oo 

P = (-% + /i) d=^ d=-z (8) 



dz dz 

For the 0+ < < oo, we choose /i = a very large 
negative number, associated with the Fock vacuum and 
obtain, for all X3, the action including the source J: 

S{J) = d^xlt/jV^j^t/jidtlj^tljidij^JijSixs)] (9) 



where V now contains a time-dependent fi{xs) that 
jumps at = from /i_ > to /i+ ^ —00. 

The generating function of the BCS wavefunctions is 



Z{J) 



(10) 



The story is depicted in the left half of Figure [l] the 
fermions travel unsuspectingly along in Euclidean time 
xs and slam like bugs onto the windshield at X3 = 0~ 
when 5{x^)Jip kills them. 

Since il) and ij) in Eq. |9[are independent Grassmann 
variables, we integrate out ijj to obtain the effective action 
for just il) to which alone J couples: 



Aid 



For the infrared limit we keep just the Jackiw-Rebbi 
zero mode^^ of the hermitian operator 



D^D(X3) = (^3 + ^i{xM-^^ + ^(2^3)), 
that obeys P/o = 

/o(a^3) = /o(0)e/o^'''(-')'*^' 
in the mode expansion of the Grassmann field: 

i'{xl,X2,Xz) = fo{X3)lp{xi,X2)- 

This kills Sind, and upon integrating /q over x^, 
SeffiJ) = f dxidx2^{id+Jfo{0))^ 



(12) 



(13) 



(14) 



(15) 



While this is indeed the action of a chiral major ana 
fermion living in the 1 — 2 plane we are not done: we 
need to show that this fermion and this action also arise 
at the edge of the same p -\- ip system. But so far we 
have no edge! It will be introduced shortly, but first a 
summary of results on the wavefunction. 

Pfaffian Wavefunction: Integrating over tjj in Eq. [Tsj 
and suppressing the constant /o (0)) we find 



ZiJ) 



exp 



1 



Aid 



J(r') 



(16) 



The two-particle wavefunction 0(ri — r2) can be writ- 
ten in terms of many related quantities: 



g^Z(J) 
dJidJ2 



1 

2id 



1 



Zl - Z2 



(17) 



and the A^-particle wavefunction is Pf( ^.i^. ). In the 
Supplementary Material we relate Z{J) and the conven- 
tional BCS wavefunction: 

\BCS) = exp Q I ^t(^)^(^ _ ^)^t(^)^^^^^ |0) (18) 

and see that (j) = —g{ri — r2). 

The Edge: To relate Z{J) in Eqn. |9]to a problem with 
the edge we rewrite S{J) in Lorentz invariant form: 



S{J) 

71 



J (fx [* (^ - /i) + J^^'] where (19) 

V'' 



J<^(a;3)(10). 



= -CTl 73 = <^3 



7^5^(20) 

(21) 
(22) 



So{J) + Si 



ind' 



(11) 



Look at the left half of Figure 1. We see our current 
description of the superconductor: translationally invari- 
ant in the xi — X2 plane, regarded as the space in which 
the pi -\-ip2 superconductor lives, and with a jump in /i at 
"time" xs = 0. In this description, the functional integral 
is saturated by one mode /o(^3), glued to the interface, 
exactly like the electron gas at a heteroj unction. 
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Extracting H{xi^ X2) from the Lorentz invariant action 
is Uke taking the row-to row transfer matrix. To derive 
the hamiltonian that governs the column-to -column dy- 
namics, we rotate the three dimensional spacetime by 
— I around the X2 axis to obtain the view shown in the 
right half of Figure [l] The points carry the same labels 
as before but the spinor undergoes a rotation: 



(23) 



Upon performing this transformation we end up with 



5(*', J) =j d?x [*' [azdi- (Tid2- (T2d3- n] 



J^(.3)(^)' 



(24) 



which describes exactly the same p -\- ip superconductor 
but in the 2 — 3 plane (with 1 ^ 3, 3 ^ — 1) with an edge 
at xs — 0. An OL(k\^k^ term may now be added without 
affecting infrared edge correlations. This is required to 
complete our identification of regions with (and without) 
fermions with \i positive (negative). 

To see that the field that J couples to is pre- 

cisely the Major ana field that arises at the edge, consider 
solving the equation for the zero mode which follows from 



Eq. 24 on dropping all Xi, X2 dependence: 



((72^3 + /i(x3))Xo = ^ Xo(^3) 



V2 V-^ 



/o(a;3)- 
(25) 

the normalizable spinor solution indeed corresponds to 
the operator -^i^^' + ^V^'^). 

We are done, for we have shown that ^( J) is at once 
the generators of electronic wavefunction in the bulk and 
of correlation functions of the Major ana field at the edge. 

For completeness, the edge Major ana field action fol- 
lows from saturating the X3 dependence of ^' as follows: 



^'(X1,X2,X3) = ^ ^. 



/0(X3)V^'(X1,X2) (26) 



Plugging this into the action S(^' ^ J) one finds, upon 
integrating the normalized function /o(^3) over x^ 

^(^', J)^ j dxidx2 [tl^'idt/j' + J/o(0)V^'] (27) 

exactly as in Eqn{T5j for the wavefunction. 

Example 2: ^He — B in D=3+l: In a simplified 
model of superfiuid ^He — Cooper pairs have spin 
1, whose projection lies perpendicular to the momenta 
^1^19 20 rpj^^ winding of this axis around the Fermi 
surface in the weak pairing phase leads to its topolog- 
ical propertie d^^ l ^-'^ i The mean-field Hamiltonian for this 



time-reversal invariant class Dili system is^ ^ is: 



H 



/i)^ka 



(28) 



[ek ■ (t]^^> 



The (i = 3 problem is just the d = 2 problem on 
steroids: A goes from being a complex number to a 
quaternion, and the spinless fermion is replaced by a two- 
component spinor. Hence the weak-pairing wavefunction 
is g^^^.{Yij) [^^^j ^^f^"^'"^^' and the many-body wavefunc- 

tion is the corresponding Pf(^) as noted in Ref. f8^. 
The Lorentz invariant action for the wavefunction is 



S 

70 



d'^x-^ 
2 

/ 

-/ 



- ji] ^ where 

_ / icre 
~ \ iecr 



/ 

-/ 



(29) 
(30) 
(31) 



Now the and 1 directions are exchanged hy R = 
exp ] , so that J now couples to ^^'^ which 

is readily verified, as before, to be the gapless edge mode 
of the rotated theory. The action for the edge theory 
obtained by saturating with the zero mode is 



/ 



rp = ^^(-a2) (32) 



Example 3: We could equally well go down a dimen- 
sion, to a spinless p-wave superconductor \n d = 1 ^ 
where A = k^^ which is also related to the quantum Ising 
model, via the Jordan- Wigner mapping. The edge the- 
ory is + 1 dimensional, corresponding to a Majorana 
zero mode, with Lagrangian C = ^i/jOxip- 

Fractionalized Topological Superconductors: We con- 
struct a fractionalized superconducting phase in D=2+l 
that bears the same relation to the p -\- ip superconduc- 
tor as the Laughlin m = 3 quantum Hall state bears 
to the integer Hall effect. Consider splitting the elec- 
tron operator at each site into three fermions ('partons') 
Cr = ifirf2rf3r and 4 ^/ir/ir/lr^^ ^ith the follow- 
ing p ^ ip mean field action for the partons: 



S{J) 
Co 



d^x[Co + iJhf2f3 



-J 

1 ^ 

= {fa fa) 



(33) 



a=l 



idi - 82 -ds - /i 



fa 
fa 



When the gauge theory is in a deconfined phase, the 
partons accurately describe the low energy dynamics. 
The SO (3) symmetry of the action, a remnant of the 
SU(3) gauge redundancy implied by = ^/ir/2r/3r^5 
is the gauge symmetry here. When the gapped bulk is 
integrated out, it generates an S0(3)i (or equivalent ly 
811(2)2)^^ Chern-Simons term which renders the gauge 
field massive thereby liberating the partons with the ac- 
tion in Eq. [33) 

Emptying out the electrons requires removing the / 
fermions, hence the strong pairing phase of the /s, where 
their chemical potential is taken to be large and negative, 
corresponds to the Fock vacuum. The electron correla- 
tors involve products of three parton correlators each in 
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a p+ip state, so the electronic wavefunction is: 





1 


jpf 





Equivalently we can start from the edge where the 
three Major ana modes are massless by gauge symmetry 
and have no relevant short range interactions in three 
spacetime dimensions. Long range gauge interactions do 
not exist due to the Chern-Simons term. Consequently 
the bulk effective action must also be described by three 
noninteracting fermions. 

Summary: We have explained why the electronic 
wavefunctions in the bulk coincided with the massless 
Majorana correlation functions at the edge in certain 
problems. We first wrote Z(J) = {$\e^'^\BCS) as a 
path integral in which the chemical potential abruptly 
jumped at in Euclidean time. Dropping the '/c^' terms 
which determined boundary conditions on /i, we obtained 
a Lorentz invariant action. Upon rotation by 7r/2 the 
same action described a system that had an edge and 
Z{J) had meanwhile morphed into the generating func- 
tion for edge correlations. In general, rotating axes will 
relate bulk wavefunctions to the edge correlations of a 
different (possibly unnatural) problem. The examples 
considered here are self-dual in this respect. 

Our analysis holds in many dimensions and applies 



to fractionalized cases as well, as long as varying ji can 
change the topology. This is possible in the Altland- 
Zirnbauer classificatioiP for models in class D in d=l and 
d=2 (like p+ip), in class C in d=2 (like d+id) and class 
Dili in d=2,3 (He-3 B phase) but not for classes like CI 
in d=3,'^ which additionally rely on band topology of the 
weak pairing Fermi surface. We are currently modifying 
our derivation for Laughlin quantum Hall states, where 
II couples to a conserved charge. 

The entanglement spectrum of the bulk seems to de- 
termine the edge theory^^'^^', which we now relate back 
to the bulk wavefunction. Since the entanglement of a 
gapped phase appears from near the cut, the entire hulk 
wavefunction must he coded holographically in every d—1 
dimensional sliver prohed in the entanglement analysis. 

Previously, the connection between edge states and 
bulk wavefunctions h as pla yed an important role identi- 
fying new FQH states^^^. Our work suggests a similar 
approach could be fruitful in identifying interacting topo- 
logical phases in D=3+l. 
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I. SUPPLEMENTARY MATERIAL 

Suppose we are given a Major ana hamiltonian 



wherein 



{*i,*i}=<5i. 



(35) 



(36) 



where z,j subsume all labels, spatial and internal. If the 
labels are continuous, the Dirac 5 should be used and 
derivatives d viewed as antisymmetric matrices. 
By definition the Grassmann integrals are 

j i/jdi/j = 1 J l-dtl; = (37) 

The Euclidean path integral corresponding to h is 



(38) 



Since coherent states do not exist for Majorana oper- 
ators (which square to ^ and not 0), one way to derive 
this result is to first form Dirac operators from pairs of 
Majorana operators, use fermion coherent states for the 
former to obtain a path integral and then undo the trans- 
formation back to Majorana fields. 

The Gaussian integral that is repeatedly used is 



Z{J) = |ei 



[dx] = Pf{A)e 



(39) 



where J and x ^i"^ 2 A/" -component Grassmann vectors. 
The two-point correlator is 



d^Z{J) 



dJadJb 



(40) 



J=0 



The minus sign in the last term can be avoided if the 
exponent is written as e~2X^x 

Higher correlators are given by Pfaffians. 



where, in the last step we have resorted to a compact no- 
tation and inserted the following resolution of the identity 
in terms of Grassmann coherent states: 



(42) 



and where it is understood for example that 



Yl \i^{x)) [##] = Yl [#(^)#(^)] (43) 



It is important to remember that ip and ip are indepen- 
dent and dummy variables. Using the defining property 
of coherent states 



(44) 



in Eq. 41 we find 



Z{J) = / [d^d^] e-'^'/'e^V'es'^s'? 



(45) 



where we have used the fact that 



since at each site 



= \o)-m 



(46) 



(47) 



and |0) = |0) ® |0) €5 ....|0). Doing the integrals over ^ 
and xp, we find 



Z{J) = e^sJ^ 



(48) 



A. PfafRan wavefunctions 



The pair wavefunction is 



Let us put these ideas to work in deriving the many- 
body wave functions from the second quantized BCS 
state. 

Consider the generating function of wavefunctions for 
any number of particles from which the wavefunctions 
can be obtained by differentiating with respect to the 
Grassmann source J{x) 



Z{J) = 



10) (41) 



(t){xi,X2) 



dJ{xi)dJ{x2) 



-g{xi-X2) (49) 



J=0 



Higher correlations follow from Wick's theorem. For 
example 

^(xi,X2,X3,X4) = g{xi-X2)g{xs-X4) 

- g{xi - X3)g{x2 - x^) 

+ g{xi - x^)g{x2 - X3). (50) 



